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In this paper, we first calculate the realistic Earth matter effects in the detection of type II super- 
nova neutrinos at the Daya Bay reactor neutrino experiment which is currently under construction. 
It is found that the Earth matter effects depend on the neutrino incident angle 0, the neutrino mass 
hierarchy Am 2 ,!, the crossing probability at the high resonance region inside the supernova, Ph, 
the neutrino temperature, T a , and the pinching parameter in the neutrino spectrum, n a . We also 
take into account the collective effects due to neutrino-neutrino interactions inside the supernova. 
With the expression for the dependence of Ph on the neutrino mixing angle #13, we obtain the 
relations between #13 and the event numbers for various reaction channels of supernova neutrinos. 
Using these relations, we propose a possible method to acquire information about #13 smaller than 
1.5°. Such a sensitivity cannot yet be achieved by the Daya Bay reactor neutrino experiment which 
has a sensitivity of the order of #13 ~ 3°. Furthermore, we apply this method to other neutrino 
experiments, i.e. Super-K, SNO, KamLAND, LVD, MinBooNE, Borexino, and Double-Chooz. We 
also study the energy spectra of the differential event numbers, dN/dE. 

PACS numbers: 14.60.Pq, 13.15.+g, 25.30.Pt, 26.30.-k 



I. INTRODUCTION 

Observable effects of supernova (SN) neutrinos in underground detectors have been a subject of intense investigation 
in astroparticle physics, both on general grounds and in relation to the SN event like 1987A In particular, flavor 
oscillation in SN may shed light on the problem of neutrino masses and mixing by means of the associated matter 
effects. Several neutrino laboratories, including the Daya Bay reactor neutrino underground laboratory [2j which is 
under construction, can be used to detect possible neutrino events from an SN explosion and serve as the SN Earth 
Warning System Q . Hence theoretical prediction for the detection of SN neutrinos in the Daya Bay and other neutrino 
experiments is very desirable. 

In Ref. the authors studied the matter effects, including both the Earth and SN, on the detection of neutrinos 
from a type II SN explosion. The results of the study were applied to the Daya Bay reactor neutrino experiment which 
measures the mixing angle #13 down to sin 2 2^13 = 0.01. In that paper, a simplified picture for the matter density 
of the Earth was used, i.e. the mantle-core-mantle picture in which p = Ylgjcrr? (p is the Earth matter density) for 
the core and p = 5g/cm 3 for the mantle, where the core radius and the thickness of the mantle are half of the Earth 
radius, respectively. The energy spectra of neutrinos from the cooling stage were assumed as exactly thermal and the 
average temperature for the neutrinos, T Ve = 3.5MeV, T Vc = bMeV, and T Vx = 8MeV, were used [y x denotes v^, v^, 
v T , or D T ). Due to the large mixing angle solution of neutrinos [f|, the crossing probability at the low resonance region 
inside the SN, Pl, is zero. Since the crossing probability at the high resonance region, Ph, depends on the unknown 
neutrino mixing angle $13, two extreme cases, i.e. Ph being and 1, were studied. In the present work, we improve 
the investigation in Ref. [J] by generalizing the simplified mantle-core-mantle picture for the Earth matter density to 
the realistic density profile including the effects due to variations of the neutrino temperature T a and the pinching 
parameter r\ a in the neutrino energy spectra. We also include the collective effects arising from neutrino-neutrino 
interactions in the SN. Furthermore, the target material is taken to be the Linear Alkyl Benzene (LAB), which has 
been selected as the main part of the liquid scintillator in the Daya Bay experiment instead of C§H\i as used in 
Ref. [|. 

In the realistic case, the matter density of the Earth changes with the depth continuously [6(, the energy spectra 
of neutrinos from the cooling stage are not exactly thermal and has a non-vanishing pinching parameter, and the 
neutrino temperature may vary in some ranges [f|[7]][|j]. Furthermore, since neutrino and anti- neutrino densities are 
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very high near the neutrino-sphere, the neutrino-neutrino interarcions become significance. This leads to collective 
effects in the SN which affect the neutrino spectra @-[l6|]. One main purpose of the present work is to investigate the 
effects of all these factors on the detection of SN neutrinos. In order to compare with Ref. [ij], we will first calculate 
the effects of the realistic distribution of the Earth matter density on the detection of SN neutrinos in the simple case 
T„ e = 3.5-MeV, T Pe = 5MeV, T x = 8MeV, in the absence of pinching parameters, i.e., r\ Vc = = r\ Vx = 0. Then 
we study the effects due to the variation of neutrino temperatures and the inclusion of the pinching parameters in the 
energy spectra. In the future, if the values of T a and r\ a can be constrained more precisely, the event number of SN 
neutrinos can be more accurately predicted. 

The other main purpose of the present paper is to propose a possible method to acquire information about $13 below 
the sensitivity of the Daya Bay experiment, i.e. 3°, through the detection of SN neutrinos. This can be achieved 
since at the high resonance region inside the SN, the crossing probability Pji is very sensitive to #13. Analytical 
expressions for transition probabilities between different neutrino flavors have been investigated in the literatures. 
Using the Landau- Zener formula [Tt]], the crossing probability, Pq, which is the neutrino jumping probability from 
one mass eigenstate to another at the resonance region, were also calculated for specific matter density distributions, 
including the linear [l8| [2(| [2l[ , exponential [221 ]. hyperbolic tangent [23|, and 1/r (24[[25[ density distributions. 
Using a similar method, the crossing probabilities at the high and low resonance regions inside the SN, Ph and Pl, 
were given in [25] 26] [27]. The neutrino event numbers N which we will obtain depend on Pjj. Using the relation 
between Ph and #13, we can predict the SN neutrino event numbers TV as a function of #13. With this, we give a 
possible method to acquire information about # 13 smaller than 1.5° by measuring the ratios of the event numbers 
of different flavor of SN neutrinos, e.g., the ratio of the event numbers of v e and v e . This method is applied to the 
existing neutrino experiments and experiments under construction, including Daya Bay, Super-K, SNO, KamLAND, 
LVD, MinBooNE, Borcxino, and Double-Chooz. 

Finally, the energy spectra of the differential neutrino event numbers, dN/dE, where E is the neutrino energy, will 
be studied. We will make predictions for dN/dE as a function of E for different neutrino Earth incident angles 9 and 
for different values of #13. Then we discuss some interesting properties of these predictions. 

The paper is organized as follows. In Sec. II, we review the necessary formulas for the detection of SN neutrinos 
on the Earth. In Sec. Ill, we apply these formulas to the Daya Bay experiment, where the realistic Earth matter 
density, the collective effects, the variation of the neutrino temperature, and the pinching of the neutrino spectra are 
considered. In Sec. IV, the relation between the neutrino event number N and #13 is discussed and a method to 
obtain information of about #13 smaller than 1.5° is proposed. This method is applied to the Daya Bay experiment 
as an explicit example. In Sec. V, the method is applied to several existing main neutrino experiments, including 
Supcr-K, SNO, KamLAND, LVD, MinBooNE, Borcxino, and another experiment under construction Double-Chooz. 
In Sec. VI, we study the energy spectra of the differential neutrino event numbers (dN/dE) and discuss some of their 
properties. Finally, a summary is given in the concluding section VII. 



II. FORMULAS FOR DETECTION OF SN NEUTRINOS ON THE EARTH 



The SN explosion is one of the most spectacular cosmic events and a source of new physical ideas. A broad area 
of topics of fundamental physics can be studied by the observation of SN. In the core collapse of SN, a vast amount 
of neutrinos are produced in two bursts. In the first burst which lasts for only a few milliseconds, electron neutrinos 
are generated via the inverse beta-decay process which leads to a neutron rich star. In the second burst which lasts 
longer, neutrinos of all flavors are produced via the nucleon-nucleon bremsstrahlung process, v c v c annihilation, and 
e+e- annihilation [28j [29t [30j . 

When the SN neutrinos of each flavor are produced they are approximately the effective mass cigenstates due to 
the extremely high matter density environment. While they propagate outward to the surface of the SN they could 
experience collective effects arising from neutrino- neutrino interactions [9|-flH| and the well-known Mikheyev-Smirnov- 
Wolfenstein (MSW) effects [l|j]-[25|]. After travelling the cosmic distance to reach the Earth, the arriving neutrinos 
are mass eigenstates, which then oscillate in flavors while going through the Earth matter. Therefore, we also have 
to consider the Earth matter effects [4] [5] 31] 32] [33] when we compute the event numbers of the various flavors of 



neutrions. 

In the Day Bay experiment, SN neutrinos undergo the following reactions in the detector. 
(1) v e + p reaction 

The large cross section, low threshold, and abundance of target protons make this the dominant channel for the 
detection of SN neutrinos. The inverse beta-decay process, 
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has a reaction threshold [34| [35j 

E th = 1.80MeV. (1) 
At low energies we approximate the large cross section as [35[ [36| 

a{D e p) = 9.5 x 10~ 44 (E{MeV) - 1.29) 2 cm 2 . (2) 

(2) v — e~ scattering 

The neutrino-electron scattering, 

v e + eT -> v e + e~ (CC and NC), 
D e + e~^D e + e~ (CC and NC), 
v a + e~ — > z/ Q + e~ (NC), 
£"a + e~ — * v a + e~ (NC), a = fx,T 

where CC and NC stand respectively for the changed current and neutral current interactions, produce recoil electrons 
with energy from zero up to the kinematics maximum. In our rate calculation we integrate over the ranges of the 
electron recoil energies. Then, the total cross sections for the neutrino-electron scattering are linearly proportional to 
the neutrino energy, and have the following forms [35l | : 

u[y e c -> v e e) = 9.20 x 10" 45 E{MeV)cm 2 , 
<j(D e e -> v e e) = 3.83 x 10~ 45 E(MeV) cm 2 , 
a{v^ r c -> v^ T e) = 1.57 x 10" 45 i;(A/eV)cTO 2 , 

<7(^ )T e -> P MjT e) = 1.29 x 10" 45 £:(7\/eV)cTO 2 . (3) 

(3) 12 C reactions 

For the neutrinos and 12 C system, there are two charged-current and six neutral-current reactions. The effective 
cross sections are obtained by scaling the experimentally measured energy values from the decay of the muon at rest 
to the energy scale for SN neutrinos. In this way, one can obtain the following effective cross sections 35] [36] 37]: 

Charged-current capture of v e : 



D e + 12 C^ 12 B + e + , E th = U.39MeV, 
12 B -^ 12 C + e- + D e , 

(a( 12 C(v e , e+) 12 B)) = 1.87 x 10 - 42 cm 2 . (4) 



Charged-current capture of v e : 



v e + 12 C -> 12 N + e" , E th = 17.teMeV, 



12 N -^ 12 C + e+ + j/ e , 

(a( 12 C(iy e , e~) 12 N)) = 1.85 x 10" 43 cm 2 . (5) 
Neutral-current inelastic scattering of v a or v a where a = e, /x, r: 

z/„+ 12 C^ 12 C+v' a , E th = 15.11MeV, 



9 a + 12 C^ 12 C*+P„, Bo, = 15.11MeV, 



i2 C * ^12 c , + 7j 

((t(i/ 12 C)) = 1.33 x 10" 43 cm 2 , 
(<r{vl 2 C)) = 6.88 x 10~ 43 cm 2 , 

H^(^) 12 C)) = 3.73 x 10~ 42 cm 2 , x = fx, r. (6) 

The effective cross sections in Eqs. (J3J-([6j are given for SN neutrinos without oscillations [3][13- When neutrino 
oscillations are taken into account, the oscillations of higher energy v x into v e result in an increased event rate since the 
expected v e energies are just at or below the charged-current reaction threshold. This leads to an increase by a factor 
of 35 for the cross section (a( 12 C{v e , e~) 12 N)) if we average it over a v e distribution with T = 8MeV rather than 
3.5MeV. Similarly, the cross section (a( 12 C(D e , e + ) 12 B)} is increased by a factor of 5. For the case of neutral-current 
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inelastic scattering of v a or v a , when the oscillations of higher energy v x into v e and v x into 9 e are taken into account, 
the cross section (a(vl 2 C)) is increased by a factor of 28 and the cross section (a{v\ 2 C)) is increased by a factor of 5. 

There will be several detectors located at the near and far sites at the Daya Bay experiment. In Ref. [H, the 
authors take the liquid scintillator to be mainly CgHi2 and the total detector mass to be 300 tons. The Daya Bay 
Collaboration has decided to use LAB as the main part of the liquid scintillator and the total detector mass is about 
300 tons. LAB, which has a chemical composition including C and H, is a mixture of the mono-alkyl benzene with 9 
to 14 carbon atoms in the side chain. The main LAB components contain 10 to 13 carbon atoms in the side chain. 
Therefore, LAB can be approximately expressed as CqH§ — C n i?2n+i where n = 9 <~ 14. In our calculation, the ratio 
of the numbers of C and H, Nc/Nh, is about 0.6. Then the total numbers of target protons, electrons, and 12 C are 

= 2.20 x 10 31 , 
A^ e) = 1.01 x 10 32 , 

N^ C) = 1.32 x 10 31 . (7) 

Neutrinos from an SN may travel through a significant portion of the Earth before reaching the detector and are 
therefore subject to the Earth matter effects. Suppose a neutrino reaches the detector with the incident angle 9 as 
indicated in Fig. 1, then the distance the neutrino travelling through the Earth is 



L= (-R + h)cos6+yR 2 -(R-h) 2 sin 2 6, (8) 

where h is the underground depth of the detector and R (6400/cm) is the radius of the Earth. At the Daya Bay 
experiment, h w 0.4fcm. Let x be the distance that the neutrino travels into the Earth, then the distance of the 
neutrino to the center of the Earth, x, is given by 



x = yJ{-R + h) 2 + (L — x) 2 + 2(R -h)(L- x) cos 9. (9) 

In the following, we will calculate the event numbers of SN neutrinos that can be observed through various reaction 
channels "i" at the Daya Bay experiment. This will be done by integrating over the neutrino energy E, the product 
of the target number Nt, the cross section of each channel a, and the neutrino flux function F® (E)/AitD 2 , 

N a (i) = N T J dE ■ a(i) ■ ■ FS, (10) 

where a stands for the neutrino or antineutrino of a given flavor, D is the distance between the SN and the Earth, 
and the index i represents different channels through which SN neutrinos are observed. 

For the neutrino of flavor a, the time-integrated neutrino energy spectra can be described by the Fermi-Dirac 
distribution (we consider the case where the spectra of neutrinos from the cooling stage are not exactly thermal) 

F (°)(E) = ^ - (11) 

where T a is the typical temperature of the neutrino 01 1 

T Vc = 3 - AMeV, T Pc = 5 - 6MeV, T Vir = 7- 9MeV, {v x = v^v T ,v^v T \ (12) 

and r\ a is the pinching parameter of the spectra (rj a > 0) to represent the deviation from being exactly thermal. The 
values of rj a for v x and v x [x = /i, t) are the same (which will be denoted as t] Utc in the following) since they have the 
same interactions, and are in general different from r\ Vii or r)p r . The values of rj a need not be constant throughout the 
cooling stage, and are typically 3 

T]v c tt 3-5, « 2.0 - 2.5, ^.RiO-2, {v x = v T , 9^, D T ). (13) 

In Eq. PT|) . F a j, where j is an integer, is defined by 



rJ 



Faj = / 7 f — rdx, (14) 

3 J exp{x-i la ) + l 
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and Na^ is the total number of the neutrinos of flavor a, 

where the average neutrino energy is (Ea* 1 ) = T^-T a and the luminosity La^ is related to the total energy release 
during the SN explosion, E s °^, through the following equation: 

4°> = ^41 (i6) 



In the numerical calculations below we take 



E { °1 = 1.97 x 10 59 MeV, (17) 



and the distance D to be 10 kpc =3.09 x 10 22 cm [3^{|4Q{. Using the above formulas, the energy spectrum function 
can be rewritten as 

T W p2 

4 0) (^)= F V , F/T T-T- (18) 

F a3 T4 exp (E/T a - r] a ) + 1 

In order to obtain the neutrino energy spectrum function at the detector, the collective effects [3]-[lIj], the MSW 
effects [lI|]-[iH: and the Earth matter effects dim (HEHI should be considered. Let P H (P L ) be the crossing 



probability at the high (low) resonance regions inside the SN, P„„ represent the collective effects of neutrino-neutrino 
interactions which is a stepwise flavor conversion probability of neutrino at a critical energy Eq, and Pi e (i = 1,2,3) 
be the probability that a neutrino mass eigenstate vi enters the surface of the Earth and arrives at the detector as an 
electron neutrino v e . Then the flux of v e at the detector, denoted as , can be written as 

i 

where Fi is the flux of Vi at the Earth surface, in either the normal or inverted hierarchy. Pj e is the probability of the 
ith mass eigenstate contained in v e and obeys the unitary condition P ie = 1. 

A significant amount of studies on the collective effects of neutrino- neutrino interactions, including simulations, 
have been made by a number of authors, e.g ., Duan et al. HI [13] 0] , Dusgupta et al. [13], Raffelt et al. [l3j[l^]. 
Esteban-Pretel et al. [13] > and Fogli et al. [l6j ] . In order to obtain a simple expression of the stepwise flavor conversion 
probabilities P vv for the neutrino and P vv for the antineutrino, we take a constant matter density and box-spectra 
for both the neutrino and antineutrino [3] ■ An analysis of the collective effects in the case of three flavors has been 
made in [13] , and it allows us to characterize the collective oscillation effects and to write down the flavor spectra of 
the neutrino and antineutrino arriving at the Earth. Following [l2| . we have P vv = P vv = 1 in the case of normal 
hierarchy; and P vv = 1, while 

in the case of inverted hierarchy, where Ec = 7MeV [l2| [f3] ■ 

Due to the large mixing angle solution of the neutrino mixing, the crossing probability at the low resonance region 
inside the SN vanishes, Pl = Pl = 0. Also, for very small sin $13, we can neglect the contributions from P^ e and 
P'ie [H|. Therefore, after a straightforward calculation the following results for the fluxes at the detector can be 
obtained: 

F£W = P 2e P H F^ + {l-P 2e P H )F^, 

<w = {1 -P 2e)F w + p 2eF (^ 

2F?W = (l-P 2e P H )F^ + (l + P 2e P H )F£\ 

2P£ W = P 2e F^ + (2-P 2e )F^\ (21) 
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2F, 



D(I) 



for the normal hierarchy (Atoj^ > 0), and 

'P2eF^> + (1 - P 2 e)F^ , (E < E C ) 

F%\ (E>E C ) 
F° (I) = P H (l-P 2e )F^ + (l + P 2e P H -P H )F^, 

' (1 - P2e)F^ + (1 + P 2 e)F^\ (E < E C ) 

F^+f£\ (E>E C ) 
2F° (I) = (1 + P 2e P H - P H )F^ + (1 + P H - P2eP H )F^ , (22) 

for the inverted hierarchy (Am^ < 0). 

Let us remark the important result that an unit flavor conversion probability means the absence of the collective 
effects. Hence as indicated in Eq. (f2"Tj) the final SN neutrino and antineutrino fluxes in the normal hierarchy are not 
modified by the collective effects. Neither are the antineutrino fluxes in the inverted hierarchy Only the neutrino 
fluxes in the inverted hierarchy and for E > Ec arc modified as indicated in Eq. (|22|) . 

The probability Pi e (i = 1,2,3) has been calculated in Ref. [3l|, in particular, 

P 2e = sin 2 6 12 + \ sin 2 29 12 f * dxV(x) sin 0™ , (23) 

2 JXQ 

where 6*12 = 32.5° [41j . V(x) is the potential that the neutrino experiences in the Earth, and * s defined as 



dzA m (x), (24) 



where 

i 

2E 

In Eq. ([25]), Am^ = 7.1 x lO~ 5 eV 2 and e(x) is defined as 



A m (x) = \j (cos 29 12 - e{x)Y + sin 2 29 12 . (25) 



, , 2EV(x) 

For a typical neutrino energy E = lOMeV, e is less than 0.13 Therefore, we neglect contributions of 0(e 2 ). 

In the Earth the potential V(x) is ^/2GfN,,(x) for an electron neutrino and —^GpN e (x) for an electron anti- 
neutrino where Gf is the Fermi constant and N e (x) is the electron number density in the Earth matter. Let the 
matter mass density inside the Earth be p{x). For nuclei of equal number of protons and neutrons the electron number 
density is 

N e (x) = p(x)/(m p + m n ), (27) 

where m p and rrin&ie respectively the proton and neutron masses. The realistic matter density inside the Earth is 
shown in Table I [6| . 

Since the crossin g pr obability at the high resonance region Pjj depends on the neutrino mixing-angle #13, which 
is unknown [28] 41] [42], we will first consider two extreme cases, Ph = (for the pure adiabatic conversion) and 
Ph = 1 (corresponding to a strong violation of the adiabatic condition). This enables us to estimate the range of the 
total neutrino event numbers to be given in Sec. III. Then in Sec. IV, we will obtain the expression of the crossing 
probability Ph and work out the relation between the event number N and #13. From appropriate ratios of event 
number which change with # 13 , we will be able to obtain some information about # 13 smaller than 1.5° by detecting 
SN neutrinos. 



III. NUMERICAL RESULTS FOR EARTH MATTER EFFECTS IN THE DETECTION OF SN 

NEUTRINOS 



In this section, we calculate the realistic Earth matter effects in the detection of type II SN neutrinos at the Daya 
Bay experiment and give the numerical results for the two extreme cases of Ph = 1 and 0. A summary of the results 
is given in Table II. 



7 



A. A simple case T Vr = 3.5MeV, T„ c = 5MeV, T x = 8MeV, T] Ve = r)v e — rj^ = 
We consider, similar to Refs. 

[lH, the sim P le 

case of T Ve = 3.5MeV, T Pe = 5MeV, T Vx = 8MeV (v x denotes v^, 
u„, v T , or D T ) and the neutrino spectra from the cooling stage to be exactly thermal, i.e., r) Vll = r]p e = r\ Vx = 0. The 
realistic Earth matter density we use is given in Table I Q . The numerical results for the event numbers for various 
reaction channels are given in Fig. 2. 

The inverse beta-decay v e + p — > e + + n has the largest Earth matter effect among all the channels. It can be seen 
from Fig. 2(a) that the maximum Earth matter effect appears at around 9 ~ 94° and the effect is about 6.82% in the 
cases of Ph = 0, 1 (normal hierarchy) and Ph = 1 (inverted hierarchy). When the neutrino incident angle becomes 
larger than about 103°, the event number increases very slowly. In the case Ph = (inverted hierarchy), the event 
number is independent of the incident angle. 

For the neutrino-electron elastic scattering, we can see from Fig. 2(b) that the total event number is much smaller 
than that of the inverse beta-decay. The maximum Earth matter effect appears at 9 ~ 93° and the amount is as 
large as 2.03% for Pjj = 1 (normal hierarchy), 0.36% for Ph = 1 (inverted hierarchy), 0.36% for Pjj = (normal 
hierarchy), while there is no Earth matter in the case Pjj = (inverted hierarchy). When the incident angle becomes 
larger than about 100°, the total event numbers increase very slowly for all the above four cases. The smallness of 
the event number of this channel is due to the small cross sections of the neutrino-electron elastic scattering given in 
Eq. ©. 

For the neutrino-carbon scattering, it can be seen from Fig. 2(c) that the maximum Earth matter effect appears at 
9 ~ 93° and the amount is as large as 1.97% for Ph = 1 (normal hierarchy), 0.99% for Ph = 1 (inverted hierarchy), 
0.99% for Ph = (normal hierarchy), while there is no Earth matter in the case Ph = (inverted hierarchy). When 
the incident angle becomes larger than about 100°, the total event numbers increase slowly for all the four cases. 

B. Realistic Earth matter effects for different T a and r\ a 

In this subsection, we consider the realistic Earth matter effects due to the variations of both the neutrino tem- 
peratures T a and the pinching parameters r\ a . In the following, we will calculate the SN neutrino event numbers for 
the limiting walue of T a and rj a given in Eqs. (|12p and Q13p . In other words, we consider the event number of SN 
neutrinos in the following two extreme cases: 

• T„ e = AMeV, T Pe = 6MeV, T Vx = 9MeV, r] Ve = 5, 7 ?Pe = 2.5, r] Vx = 2; 

• T v . = 3MeV, T Pe = 5MeV, T Vx = 7MeV, r) Ue = 3, jjp. = 2, rj Ux = 0. 
The results are shown in Figs. 3-4. 

(1) The case T„ e = 4MeV, = 6MeV, T Vx = 9MeV, r) Ve = 5, tj^ = 2.5, t] Vm = 2 

For the inverse beta-decay v e +p — > e + + n, it can be seen from Fig. 3(a) that the maximum Earth matter effect 
appears at around 9 ~ 95° and the effect is about 6.98% in the cases where Ph = 0,1 (normal hierarchy) and Ph = 1 
(inverted hierarchy). When the incident angle becomes larger than about 103°, the event number increases very 
slowly. In the case Ph = (inverted hierarchy), the event number is independent of the incident angle. 

We can see from Fig. 3(b) that the total event number of the neutrino-electron elastic scattering is much smaller 
than that of the inverse beta-decay. The maximum Earth matter effect appears at 9 ~ 94° and the amount is as 
large as 2.48% for Ph = 1 (normal hierarchy), 0.43% for Ph = 1 (inverted hierarchy), 0.43% for Ph = (normal 
hierarchy), while there is no Earth matter in the case Ph = (inverted hierarchy). When the incident angle becomes 
larger than about 100°, the total event numbers increase very slowly for all the above four cases. 

For the neutrino-carbon scattering, it can be seen from Fig. 3(c) that the maximum Earth matter effect appears at 
9 ~ 93° and the amount is as large as 2.86% for Ph = 1 (normal hierarchy), 1.43% for Ph = 1 (inverted hierarchy), 
1.43% for Ph = (normal hierarchy), while there is no Earth matter in the case Ph = (inverted hierarchy). When 
the incident angle becomes larger than about 100°, the total event numbers increase slowly for all the four cases. 

(2) The case T Ve = 3MeV, = 5MeV, T Vx = 7MeV, ^ = 3, ^ = 2, r] Ux = 

Similar to the above case, for the inverse beta-decay v e + p — > e + + n, it can be seen from Fig. 4(a) that the 
maximum Earth matter effect appears at around 9 ~ 94° and the effect is about 3.69% in the cases where Ph = 0, 1 
(normal hierarchy) and Ph = 1 (inverted hierarchy). When the incident angle becomes larger than about 103°, the 
event number increases very slowly. In the case Ph = (inverted hierarchy), the event number is independent of the 
incident angle. 

For the neutrino-electron elastic scattering, which is plotted in Fig. 4(b), the maximum Earth matter effect appears 
at 9 ~ 93° and the amount is as large as 1.59% for Ph = 1 (normal hierarchy), 0.21% for Ph = 1 (inverted hierarchy), 
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0.21% for Ph = (normal hierarchy), while there is no Earth matter in the case Ph = (inverted hierarchy). When 
the incident angle becomes larger than about 100°, the total event numbers increase very slowly for all the above four 
cases. 

Also, for the neutrino-carbon scattering, it can be seen from Fig. 4(c) that the maximum Earth matter effect appears 
at 8 ~ 93° and the amount is as large as 1.62% for Ph = 1 (normal hierarchy), 8 ~ 92° and 0.81% for Pg = 1 (inverted 
hierarchy), 8 ~ 92° and 0.80% for Ph = (normal hierarchy), while there is no Earth matter in the case Ph = 
(inverted hierarchy). When the incident angle becomes larger than about 100° , the total event numbers increase 
slowly for all the four cases. 

The above results can be understood by considering the fact that the oscillation behavior is determined by the 
factor l\m\ l (eV 2 )L(m)/E(MeV), where L is the distance that the neutrino travels in the Earth as given in Eq. (jHJ). 
When 8 < 90° this distance is smaller than 10 km and hence the amount of the Earth matter effects are very small. 
When 8 increases beyond 90° this distance exceeds 100 km, the Earth matter effect becomes greater and reaches a 
maximum value for 91° ~ 95°. When 8 is more than 100°, the distance that neutrino travels in the Earth is greater 
than 2000 km, then there could be many oscillations in and hence the averaging Earth matter effect is smaller 
than the maximum value. 

In Table II, we list the incipient values of the SN neutrino event numbers where the neutrino incident angle 8 is zero 
and hence practically the vacuum event numbers, the minimum SN neutrino event numbers and hence the maximum 
Earth matter effects, the incident angles at which the maximum Earth matter effects appear, and the Earth matter 
effects for all the cases considered. It can be seen that the realistic Earth matter effects are much smaller than those 
given in Ref. Q where the mantle-core-mantle approximation for the Earth matter density was used. When both 
T a and r\ a vary in their ranges, the event numbers of SN neutrinos and the Earth matter effects will vary all the 
three reactions. More precise values of T a and r\ a will help obtain more reliable event numbers and the Earth matter 
effects. There are also some characters which can be seen from Table II: 

• In the two cases: Ph = 1 (inverted hierarchy) and Ph = (normal hierarchy), we have the same event numbers 
and the same Earth matter effects for all the three kinds of reactions; 

• In the case Ph = (inverted hierarchy), the event numbers are independent of the incident angle 8 (i.e. there is 
no Earth matter effects) for all the three kinds of reactions; 

• In the case Pjj = 1 (normal hierarchies), there are the largest Earth matter effects in all the four cases for the 
three kinds of reactions; 

• When T a and rj a increase, the maximum Earth matter effects in all the four cases increase for the three kinds of 
reactions. 



IV. A POSSIBLE METHOD TO ACQUIRE INFORMATION ABOUT 9xs SMALLER THAN 1.5° 

In the above section, we took the crossing probability at the high resonance region inside the SN, Ph, to be in the 
two extreme cases, i.e. Ph = and Ph = 1, due to the fact that Ph depends on the neutrino mixing angle #13 which 
is unknown. In this section, using the expression of Ph, wc will derive relations between the event numbers of SN 
neutrinos, N, and #13, and then pro pose a possible method to obtain information about #13 smaller than 1.5°. 

Considering the MSW effects jl8l ] , a number of authors including Bethe [lj| , Parke , Haxton [2l[ , Petcov et 
al. [lH, Notzold [23||, Kuo and Pantaleone (24|[25[ used the Landau-Zener formula jl7l| to calculate the crossing 
probability Pc at the resonance region inside a star. Several different density distributions in the case of two-flavor 
transitions were considered. In particular in [24|(25j the following result is given: 

p cxp(-f 7 F)-exphf 7 (^)] 
l-exp[-f 7 (^)] 



where [3 is the mixing angle and 



^sin^ 1 

' 2E cos2/3 \d\nN e /dr\ res ' V ' 



where N e is the electron density and F can be calculated by Landau's method. The expression of F was given in 
Table I in Ref. [24| and Table III in Ref. [25| for different density distributions. Consequently, if we consider the 
following electron density distribution 



(30) 
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where k is a constant and n is an integer, then 



7 = 



2\n\ 



| Am 2 1 
E 



sin 2 2/3 \ / cos 2/3 
cos 2/3 J I 2V2G F ^- 



(31) 



F = 2V( 

/-^ V 2m 



1/2 
m + 1 



(tan 2/3) 



2m 



(32) 



where 



1/n- 1 
2m 



(1/n-l)! 



(1/n - 1 - 2m)!(2m)! ! 



1/2 
m + 1 



(-1)' 



5^4 ' 



(sin< 



(2m- 1)!! 7T 
(2m)!! 2' 



(33) 



We note that our expression for J m in Eq. (|3"3"j) has the same final result as given in Eq. (B5) in Appendix B of Ref. 
[2~il ]. But the defining integrals are different. 

In order to obtain the expression for the crossing probability at the high resonance region inside the SN, Ph, we can 
use the similar method in the case of three-flavor transitions [2(| [13] . It is known that the following matter density 
profile for the SN is appropriate (2^4100: 



10 cm 



10 iu g/c 



(34) 



where the constant C depends on the amount of electron capture during the star collapse and its value is between 1 
and 15. This corresponds to n = —3 and k = C ■ 10 31 in Eq. ([30]) . then we have 



1 



10 10 MeV ( sin 3 26 



13 



E \ cos 2 26»i3 / V lel^ 2 



2/3 



For very small #13, Ph has a simpler expression: 



P H = exp(--jF). 

Since #13 is small, we need only to consider the m = term in Eq. Q32p . then 



F^2( 



In addition, since Am 2 x <C lAmJ-J, we can set 



lAm 2 ,! « |Am 2 



32 



Then we obtain 

P H = cxp I - 
In the 2a allowed ranges (U [H| [H] [47 



7T 
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10 10 MeV ( sin 3 26 u 



E \ cos 2 26»i3 



= 1. 



leV 3 



| C l/2 



2/3- 



|Am 2 2 | = 2.6 x x lO^ 2 , 

^13 < 10°, 



(35) 



(36) 
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and we will use 2.6 x 10~ 3 eV 2 for |Am| 2 | in the following calculation. 

From Eq. (|35[) . the dependence of Ph on 13 and E is shown explicitly. In order to obtain such simple expression 
of Pjj, we assume #13 to be very small and use the simplified structure model of SN in which p oc ?'~ 3 in the above 
calculations. Therefore, Eq. (|35[) is usually suitable for very small 613 (we have checked numerically that in the range 
of 0i3 we are working Eq. (|35[) is a good approximation), and the expression for Ph should be reevaluated if #13 is 
larger. 

In Fig. 5(a), we plot Ph as a function of #13 and it can be seen that when #13 is zero, Ph = 1; in the small range 
of 813 of 0° ~ 2° Ph varies rapidly in the range < Ph < 1; and when #13 is larger than 2°, Ph = 0. Therefore, it is 
possible to use the value of Ph to obtain the value of 613 when it is in the range 0° ~ 2°. In Fig. 5(b), we plot Ph as 
a function of E with 13 in the range 0° ~ 2° and again it shows that 613 can be measured by the energy spectrum 
of Ph if the value of 613 is in the range of 0° ~ 2°. 

Using Eqs. (fTOj) . ([T5jl . (f2"Tj) . (j2"2"|) , and ([55]) . we obtain the relation between the event number of SN neutrinos N 
and the neutrino mixing angle #13. The curves of N versus 613 in the range 0° ~ 3° for the Daya Bay experiment for 
the inverse beta-decay is given in Fig. 6. From this plot, we can see that when #13 < 1.5°, N is very sensitive to #13. 
When 813 > 1.5°, N is nearly independent of $13. Therefore, when #13 is smaller than 1.5°, one could acquire some 
information about #13 by detecting the event numbers of SN neutrinos. However, in practice, since we do not have 
the accurate values of the SN neutrino flux parameters T a and r\ a and they vary in some ranges, the event number 
of SN neutrinos which we will detect depend on these two unknown parameters strongly. It can be seen from Fig. 6 
that the uncertainties of the event number due to those of T a and r\ a are very large, and it is difficult to obtain the 
value of #13 from the event number of SN neutrinos. This difficulty also occurs in the neutrino-electron scattering 
and neutrino-carbon reactions. Therefore, we need to work with a quantity which changes with 613 but is insensitive 
to the values of T a and r\ a . 

In order to reduce uncertainties from unknown quantities such as luminosity Q, distance [40l |. temperatures T a , and 
pinching parameters r] a , ratios of measurable quantities can be used. For instance, [48| defines the ratio of high-energy 
to low-energy event numbers that are measurable in neutrino oscillation experiments: 

N evmt (30 < E < 70MeV) 

N event {5 <E< 20MeV) ' [ ' 

Also, many different quantities to obtain the information about 613 and the mass hierarchy are suggested in [26j]. In 
the following we choose a suitable reaction and then define a ratio of event numbers which is sensitive to the mixing 
angle #13 but only dependents on T a and rj a slightly. 

For the inverse beta-decay, there is only one flavor neutrino v e . The numerator and denominator or R of Eq. 
(j3"T|) have the same flavor and our calculation shows that it has a significant dependence on T a and r\ a . Hence R is 
not a suitable quantity to work with. In the case of neutrino-electron scattering, the event number of SN neutrinos 
detected at the Daya Bay experiment is very small as shown in the above section. Again this channel of reaction is 
not useful for our purpose. We now turn our attention to the channel of the neutrino-carbon reactions at the Daya 
Bay experiment. 

For the neutrino-carbon reactions, we define the quantity R to be the ratio which is the event number of u e over 
that of v e . Using Eqs. (TIT))) . (TT5|) . (|2~lj) , (|2"2"|) . and (f3"5)) . we can obtain the relation between R and the mixing angle 
#13. With this, we can obtain some information about 613 smaller than 1.5° since the uncertainties due to T a and r\ a 
are small. In Fig. 7, we plot R as a function of the mixing angle #13 when T a and rj a take their limiting values in the 
ranges given in (TT2")) and (fT3"|) for different incident angle 6. It can be seen from these plots that the uncertainties of 
R due to T a and r\ a are indeed not large. For 6*13 < 1.5°, R is very sensitive to 6*13. However, while 613 > 1.5°, R is 
nearly independent of #13. Therefore, when #13 is smaller than 1.5°, we may restrict the mixing angle #13 in a small 
range and get information about mass hierarchy by detecting the ratio of event numbers of SN neutrinos even though 
there are still some uncertainties due to the incident angle 6, the mass hierarchy Am^, and the structure coefficient 
C of the SN density function as given in Eq. (j34|) . Recently, there are discussions on methods for the determination 
of the neutrino mass hierarchy and the incident angle of the SN. For examples, in Ref. (49|, a method to identify 
the mass hierarchy at extremely small #13 through the Earth matter effects is g iven, and a method to determine the 
incident angle of the SN by the electron scattering events can be found in 50|. For the structure coefficient C, its 
value can be more precisely determined if we can obtain more exact density profile of the SN (45| . In the future, if the 
incident angle 6, the mass hierarchy Am^, and the the structure coefficient C can be determined, the uncertainties 
in the determination of 613 through the SN neutrino will be much reduced. 

At the Daya Bay experiment, the sensitivity of sin 2 29\3 will reach 0.01, i.e., to determine #13 down to about 3°. 
Therefore, if the actual value of 613 is smaller than 3°, the Daya Bay experiment can only provide an upper limit for 
#13. However, if an SN explosion takes place during the operation of Daya Bay, roughly within the cosmic distance 
considered here, it is possible to reach a much smaller value of #13 through the ratio of the event number of different 
flavor SN neutrinos in the channel of neutrino-carbon reactions as discussed above. 
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V. ACQUIRING INFORMATION OF 6 1S SMALLER THAN 1.5° FROM OTHER NEUTRINO 

EXPERIMENTS 



In this section, we will apply the above method to some other current neutrino experiments including Super-K, 
SNO, KamLAND, LVD, MinBooNE, Borexino, and Double-Chooz which is under construction. 

For Super-K, the target material is water, the total detector mass is 32000 tons, and the depth of the detector 
h = 2700 m.w.e. 1 [13] [48] 51} [52]. Then the total numbers of the targets (protons, electrons, and 16 O) are 

= 2.14 x 10 33 , 
N^ e) = 1.07 x 10 34 , 

N^ 0) = 1.07 x 10 33 . (38) 

First, we consider the inverse beta-decay. When the incident angle 8 = 30°, we use Eqs. (fTO]). (fl~8|) . (l2Tj) . (j22|) . and 
()35p to plot the event number of SN neutrinos at Super-K as a function of the mixing angle # 13 . The result is shown 
in Fig. 8. Similar to the Daya Bay experiment, it can be seen that the uncertainties of the event number due to the 
two quantities T a and rj a are very large, and hence we cannot make predictions on #13 from the event number of SN 
neutrinos. Therefore, we need to choose a suitable reaction at Super-K and find out a quantity which is sensitive to 
013 but only dependents on T a and rj a slightly. 

Although the event number of SN neutrinos in the channel of neutrino-electron scattering is very small at Daya Bay, 

it is much larger at Super-K due to the large number of target electrons Nj, . We define R as the ratio of the event 
number of v e to that of v e from SN. Using Eqs. (p~D|) . (fTS)) . (|21[) . (|2"2"|) . and (f3"S"|) . we can obtain the relation between 
R and the mixing angle #13. By measuring this ratio, we can acquire some useful information about #13 smaller than 
1.5° since the uncertainties due to T a and rj a are small. In Fig. 9, we plot R as a function of the mixing angle #13 
when T a and r\ a take their limiting values in the ranges (|12[) and (|13p for different incident angle 6. It can be seen 
from these plots that the uncertainties in R due to T a and r\ a are not large. When #13 < 1.5°, R is very sensitive to 
6*i 3 . However, while $13 > 1.5°, R is nearly independent of $13. Therefore, when 8 13 is smaller than 1.5°, #13 can be 
constrained in a small range by the event numbers of SN neutrinos. Furthermore, in spite of uncertainties due to #13, 
T a , rj a and C, the plots show that if R is smaller than about 2.35, the mass hierarchy must be normal while if R is 
larger than about 2.5, the mass hierarchy must be inverted. 

For SNO, the detector material was heavy water, the total detector mass was 1000 tons, and the depth of the 
detector h = 6000 m.w.e. 53| 54|. Therefore, we need to consider the reactions between neutrinos and deuterium. 



There are two changed-current and six neutral-current reactions. The cross sections are obtained by scaling the 
experimentally measured energy values from the decay of the muon at rest to the energy scale for SN neutrinos. In 
this way, one can obtain [36] 48] [55] . 5(| : 



Charged-current capture of v e or v e 



v e + d -> p + p + e~, E t h = lAAMeV, 
9 e +d^n + n + e + , E t h = 4.03M eV, 
(a(d(v e ,e-)pp)) = (3.35T£ 31 - 3.70) x 10 - 43 cm 2 , 
(a(d(v e ,e+)nn)) = (3.05T^ 08 - 7.82) x 10 " 43 cm 2 . 



Neutral-current inelastic scattering of v a or v a where a = e, n, r: 



v a +d^ n + p + v a , E th = 2.22MeV, 
& a + d^n + p + &' a , E th = 2.22Af eV, 
(a(d(v a ,v a )np)) = (1.63T£ 26 - 2.78) x 10 " 43 cm 2 , 
(a(d(v a ,v' a )np)) = (2.03T 2 05 - 3.76) x 10 ~ 43 cm 2 . 



1 m.w.e. refers to meter-water-equivalent and 1 m of rock is about 2.7 m of water. 
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For T„ c = 3.5MeV, T Pc = 5MeV, T„ x = 8MeV, 



(a{d(v e , 


e )pp)) = 


C dO w 1(1- 

5. Do X 11) 


-42 2 

cm , 


{a(d(v e , 


e Jnn)) = 


7 on \y 1 Pr 


-42 2 

cm 


(<r(d(u e , 


v' a )np)) = 


2.49 x 10" 


A2 cm\ 


(a(d(v e , 


K)np)) = 


5.12 x 10" 


A2 cm\ 


(a(d(u x , 


v 'x) n V)) = 


1.76 x 10" 


-41 2 

era , 


(a(d(v x , 


v'x) n v)) = 


1.40 x 10" 


-41 2 

era , 



As indicated in Eq. (|39[) . it can be seen that the average effective cross sections of v x and reactions arc different. 
Therefore, we need to distinguish the energy spectrum functions F$ and F^J appearing in Eqs. (|2~lj) and (|22[) . 

Similar to the neutrino-carbon scattering, the average effective cross sections in Eq. (j3"9")) are given for SN neutrinos 
without oscillations [56| . When neutrino oscillations are taken into account, the oscillation of higher energy v x into 
v & again results in an increased event rate since the expected v e energies are just at or below the chargcd-currcnt 
reaction threshold. This leads to an increase by a factor of 7 for the cross section (a(d(v e , e~)pp)). Similarly, the cross 
section (a(d(D e , e + )nn)) is increased by a factor of 3. For the case of neutral-current inelastic scattering of v a or 9 a , 
when the oscillation of higher energy v x into v e is taken into account, the cross section (a{d(y e , u a )np)) is increased 
by a factor of 7, while when the oscillation of higher energy v x into v e is considered, the cross section (a(d{p e , v a )np)) 
is increased by a factor of 3. 

Now, we define R as the ratio of the event number of v e to that of v e in the channel of neutrinos-deuterium 
scattering. Using Eqs. (fTU|) . IfTS]). (f2"T) . (|2"2"|) . and ((331) . wc can obtain the relation between R and the mixing angle 
#13. In Fig. 10, we plot R as a function of d\j, when T a and r/ a take their limiting values in (fT2")l and (fT"3")) for different 
incident angle 8. Similar to electron-carbon reactions, it can be seen from these plots that the uncertainties of R due 
to T a and r\ a are not large. When 6*13 < 1.5°, R is very sensitive to #13. However, while #13 > 1.5°, R is nearly 
independent of #13. Therefore, when # 13 is smaller than 1.5°, we can constrain # 13 in a small range by the detection 
of the event numbers of SN neutrinos in the channel of neutrinos-deuterium reactions at SNO. 

Similar to the Daya Bay experiment, the method to obtain information about #13 smaller than 1.5° can be applied 
to KamLAND [|37]|, LVD [58[, MinBooNE [59[, Borexino (6(j, and Double-Chooz [61| in the channel of neutrino-carbon 
reactions. We can also constrain #13 in a small range with the ratio of the event number of SN v e to that of v e in 
these experiments. Since the number of target carbons of Doublc-Chooz is much smaller than that of Daya Bay, the 
event number of SN neutrinos which can be detected is fewer and it will be more difficult to obtain information about 
013 by our method. 



VI. ENERGY SPECTRA OF THE DIFFERENTIAL EVENT NUMBERS 

In this section, we consider the energy spectra of the differential event numbers, referred generically as dN/dE, 
where E is the neutrino energy. From Eq. (|10[) . one can obtain 

dNa ^ = N T aU)—F D (40) 

Since the inverse beta-decay is the most important reaction among the three kinds of reactions in the Daya Bay 
experiment, we examine the energy spectra of this process, taking the simplest parameter set T Vi! — 3.5MeV, Tp e = 
hMeV , T Vx = 8MeV, r)„ e = 7?p e = r/„ T = as an example. 

Using Eqs. (f2"T) . (f2"2"]) , and (|4*U)) . we make a three-dimensional plot of dN /dE versus E in the Daya Bay 

experiment for different Earth incident angle 9 in Fig. 11(a) for # 13 = (inverted hierarchy). We can see that when 
the incident 8 increases the Earth matter effects become more and more obvious. The curve of dN /dE changing with 
E is very smooth when 8 is small. However, it becomes oscillatory when 8 is greater than 90°. Therefore, this shows 
that the Earth matter effects on dN /dE need to be considered in this case. 

Using Eqs. ©, JEJl, jUJ), (22]), (35}, and (gOJl, we plot dN/dE as a function of E for different 8 13 in Fig. 11(b) 
for the case of inverted hierarchy, where the structure coefficient of SN C = 3, and the Earth incident angle 8 = 30° 
are used. We can see that the energy spectrum changes with #13 in the range 0° ~ 1.5°. 
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VII. SUMMARY AND DISCUSSIONS 



In this paper, we have calculated the realistic Earth matter effects in the detection of type II SN neutrinos at the 
Daya Bay experiment under construction. It is found that the Earth matter effects depend on the neutrino incident 
angle 9, the neutrino mass hierarchy Am§ 2 , the crossing probability at the high resonance region inside the SN, Pff, 
the neutrino temperature T a , the pinching parameter in the neutrino energy spectra r\ a , and the collective effects of 
neutrino-neutrino interactions in the SN. We have given the event numbers that can be detected through the inverse 
beta-decay, the neutrino-electron scattering, and the neutrino-carbon scattering. We have studied the effects due to 
the variations of the neutrino temperature T a and the pinching parameter rj a in the neutrino energy spectra. 

Since neutrino flavor conversions inside the SN depend on the neutrino mixing angle #13, it is possible to get 
information about #13 by detecting SN neutrinos. In fact, there have been some general discussions on this possibility 
[Ifl [29| . We have made concrete calculations and given explicit event numbers by using the relation between the event 
numbers of SN neutrinos N and the neutrino mixing angle #13 under different scenarios of neutrino parameters. For 
613 smaller than 1.5°, we propose an approach to constrain the value of angle #13 in a small range and get information 
about mass hierarchy by measuring the ratio of the event numbers of different flavors of SN neutrinos. For the 
Daya Bay experiment, we choose the ratio of the event number of v e to that of v e in the channel of neutrino-carbon 
reactions. We have also applied this method to other neutrino detectors including Super-K, SNO, KamLAND, LVD, 
MinBooNE, Borexino, and Double-Chooz. For the Super-K, the suitable reaction is the neutrino-electron scattering 
due to its large number of target electrons. We have also shown that the neutrino-deuterium reactions at SNO may 
be used to acquire useful information about #13. 

We have studied the energy spectra of the differential event numbers, dN/dE. From the dependence of dN/dE 
on the neutrino energy E for different Earth incident angle 9, we have found that when the Earth incident angle 9 
increases the Earth matter effects become more and more pronounced. 

There are still some uncertainties in our work. In our calculations, we took the distance from the SN to the Earth 
to be 10 kpc, where the maximum of the progenitor population appears in the Milky Way (39j . The distribution of 
the SN progenitors as a function of the distance to the Earth can be found in Ref. [4Cl |. 

We let the parameters in the neutrino energy spectra (the temperatures and the pinching parameters) vary in some 
reasonable ranges. In fact, the simulations from the two leading groups, the Livcrmore group [H[ and the Garching 
group @ (which considered more reactions in their simulation and found different dominant neutrino production 
processes in the formation of the SN neutrino spectra and fluxes), led to parameters which agree within about 20- 
30%. However, their central values of SN parameters are different. 
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Figure Captions 

Fig. 1 Illustration of the path of the SN neutrino reaching the detector in the Earth. D is the location of the detector, 
9 is the incident angle of the neutrino, O is the center of the Earth, L is the distance the neutrino travels through the 
Earth, and x is the distance of the neutrino to the center of the Earth. 

Fig. 2 The event numbers observed at the Daya Bay experiment as a function of the incident angle 9 when T Vc = 
3.5MeV, T Pe = 5MeV, T Vx = 8MeV, and r\ Vc = r]^ = rj Vx = 0. (a) the channel D e + p — > e + + n; (b) the reactions 
v + e~ — ► v + e~; (c) the neutrino-carbon reactions. The solid curves correspond to Ph = 1 (normal hierarchy), the 
dashed curves correspond to Ph = 1 (inverted hierarchy), the dotted curves correspond to Ph = (normal hierarchy), 
and the dot-dashed curves correspond to Ph = (inverted hierarchy). For (a), the solid curve also corresponds to 
Ph = 1 (inverted hierarchies) and Ph = (normal hierarchy). 

Fig. 3 Similar to Fig. 2 but T a and i] a take their maximum values: T„ e = 4MeV, T Pe = QMeV, T Vx = 9MeV, 
f]v e = 5, r]y e = 2.5, T] Vx = 2. 

Fig. 4 Similar to Fig. 2 but T a and r\ a take their minimum values: T Ve = 3M eV, T^^ = 5MeV, T Vm = 7MeV, T] Vc = 3, 
r)ve = 2, rf Vm = 0. 

Fig. 5 The crossing probability at the high resonance region inside the SN: (a) as a function of the neutrino mixing 
angle 6>i3 and the solid curves, dotted curves, dashed curves correspond to E = HMeV, l6AleV, 25MeV, respectively; 
(b) as a function of the neutrinos energy E and the solid curves, dotted curves, dashed curves correspond to #13 = 
0.1°, 1°, 2°, respectively. 

Fig. 6 The event number observed in the channel v e + p — > e + + n at the Daya Bay experiment as a function of the 
neutrino mixing angle #13 when the incident angle 9 — 30°. The solid curves correspond to the normal hierarchy, 
and the dashed curves correspond to the inverted hierarchy, where "max" ("min") corresponds to the maximum 
(minimum) values of T a and r\ a - 

Fig. 7 The ratio of the event number of v e to that of P e , R, as a function of the mixing angle #13 in the channel of 
neutrino-carbon reactions at the Daya Bay experiment, (a) the incident angle 9 — 30°; (b) 9 = 90°; (c) 9 = 93°; (d) 
9 = 150°. The solid curves correspond to the normal hierarchy (max), the dashed curves correspond to the inverted 
hierarchy (max), the dot-dashed curves correspond to the normal hierarchy (min), the dotted curves correspond to 
the inverted hierarchy (min) , where " max" (" min" ) corresponds to the maximum (minimum) values of T a and r\ a . 

Fig. 8 Similar to Fig. 6 but for the Super-K experiment. 

Fig. 9 Similar to Fig. 7 but for the ratio of the event number of v e to that of v ei R, as a function of the mixing angle 
9i3 in the channel of neutrino-electron scattering at the Super-K experiment. 

Fig. 10 Similar to Fig. 7 for the ratio of the event number of v e to that of P e , R, as a function of the mixing angle 
#13 in the channel of neutrino-deuterium reactions at the SNO experiment. 

Fig. 11 In the case of inverted hierarchy, the differential event number dN/dE observed in the inverse beta-decay 
channel at the Daya Bay experiment as a function of the neutrino energy E. (a) for different incident angle 9 (#13 = 0); 
(b) for different 13 (9 = 30°, C = 3) when T„ c = 3.5MeV, = 5MeV, T Ux = SMeV, r) Ve = r?p e = r) Vx = 0. 
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TABLE I: The realistic matter density inside the Earth where r is the distance to the center of the Earth and 7? — r (km)/ 6371 
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TABLE II: Summary of the realistic Earth matter effects. N (7) represents the normal (inverted) hierarchy, 'I' ('0') represents 
Ph — 1 (0). The numbers in the columns 'Incipient' and 'Min' are the event numbers when the SN neutrino Earch incident 
angle is zero and is the angle in the column 'Angle', respectively. The column Angle' gives the angles at which the event 
numbers are minimum and the Earth matter effects are the strongest. The column 'Ratio' gives the percentages of the Earth 
matter effects. 
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